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ABSTRACT 


A progran to study stress generation mechanissis in silicon 
sheet growth was started at Mobil Solar on July 1982. The 
purpose of th<^ research is to define post-growth temperature 
profiles for the sheet that can minimize its stress during growth 
at high speeds, e.g., greater than 3 cm/min. The initial tasks 
described in this report concern work in progress toward the 
development of computing capabilities to (i) model 
stress- temperature relationships in st >ady-state ribbon growth, 
and (ii) provide a means to calculate realistic temperature fields 
in ribbon, given growth system component temperatures as boundary 
conditions. If it is determined that low stress configurations 
can be achieved, the modeling is to be tested experimentally by 
constructing low-stress growth systems for EFG silicon ribbon. 


•The JPL Flat Plate Solar Array Project is sponsored by the U.S. 
Department of Energy and forms part of the Solar Photovoltaic 
Conversion program to iniliate a major effort toward the 
development of flat plate solar arrays. This work was performed 
for the Jet Propulsion Laboratory, California Institute of 
Technology by agreement between NASA and DOE." 
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I. IHTRODUCTIOM 

A satisfactory nodal that can account for strassas ganaratad 
in silicon shaat grown at high speeds is not yet available, 
tfunerous attempts to account for residual stresses have been nade, 
but all of these suffer inadequacies in one area or another. The 
most prominent of these are the lack of a fully two-dimensional 
treatment of stress and temperature fields in the sheet width 
dimension, and the omission of plastic deformation effects. Two 
significant factors preventing the development of such a model 
have been the absence of adequate information on creep behavior of 
silicon at high temperatures that could be applicable to the sheet 
growth situation, and of experimental data on temperature fields 
and deformation modes of growing sheet that could guide the 
modeling effort. 

This report describes the work in progress under this 
subcontract to attempt to develop and test in the laboratory a 
stress- temperature field model for EFG silicon ribbon that deals 
with the above deficiencies. In one subtask, a computer code to 
predict stress- temperature field relationships in steady-state 
sheet growth is being developed at Harvard University. The stress 
state will be parameterized by a two-dimensional temperature field 
and growth speed. Incorporation of time dependent stress 
relaxation effects through a creep law is also planned to model 
the impact of plastic flow on the sheet residual stress state. A 
second aspect of the program will deal with the development of a 
model to predict the temperature field in a moving sheet from 
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givv^n tysten conponent tftBp«raturts th« ahaat anvlronnant ) , 

and study axparimental ueana to varlfy the nodal. 

The goal of the program is to combine the results of these 
two areas of study to arrive at a model that can predict 
stress- temperature field relationships in steady-state silicon 
sheet growth under realistic conditions. Minimum stress 
configurations will be sought, and an attempt to construct an EFG 
silicon ribbon growth system that can verify this model will be 
made if it appears such configurations can be achieved 
experimentally. To aid the effort to correlate modeling results 
to experimental conditions, experimental means to examine 
temperature fields and stress generating process during growth are 
also under study. 

II. PROGRESS REPOR T 

The first section that follows describes the strategy that 
has been adopted to develop a computer code to model the 
stress- temperature field relationship in a silicon sheet growing 
under steady- state conditions. This work is being done at Harvard 
University. The second section focuses on the development of a 
computer program to calculate ribbon temperatures for given 
temperature boundary conditions in an enclosure surrounding the 
ribbon, and to verify the calculation procedure experimentally. 
This part of the work is being done at Mobil Solar. 


2 




OF POOH QUAvtrr 


A. A flalyg la oi SteadyrState Ribbon Growth (J.U. 

Hutchinson and J.C. Lanbropoulos , Harvard Univarsity) 

Work in the period July 9, 1982, to October 1, 1982, was 

largely directed towards developing a computer code, for analyzing 

the stress distribution in a ribbon pulled from the melt through 

an arbitrary temperature field under steady>state conditions. The 

formulation of the problem and the numerical method for solving it 

are described below. Similar methods have been applied with 

success to analyze steady- state crack growth problems in 

elastic- plastic solids^ and in elastic-creeping solids. 

The ribbon emerges from the melt at x = 0 and is pulled in 

the positive x-direction with uniform velocity V. It has width 2h 

extending from y - -h to y = h. The temperature distribution in. 

the ribbon is denoted by T(x,y) and is assumed to be independent 

of time. 

Plane stress conditions are assumed with non- zero stresses 

<r , <7 and o' . Steady-state conditions are assumed to exist so 
XX yy xy 

that the stresses, strain- rates and velocities are independent of 

time at any spatial point (but not, of course, for a material 

element). The additional velocity components are denoted by 

(v ,v ) so that the total velocity components are (V ♦ 

X y * y 

where the additional velocities are very small compared to the 
pull velocity V. The strain-rates are related to the velocities 
by 








( 1 ) 
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Thft nat«rial is taken to be an elaetlc-ereeping aolid. The 
nuaerlcal nethod la not tied to any particular form of the 
constitutive lav. However, the law which looks nost pronlalng at 
the noment to represent existing creep and relaxation data for 
silicon has the form of the following prototype relation: 

( 2 ) 


0 . - li,.,. 4„+|bo ”*V,+ ai«, 


ij ® ■'ij" ® 

Here, the dot denotes the time-rate of change with respect to 
a material element and, for steady-state pulling through a 
time- Independent temperature distribution. 


T»VT,^ (3) 

The Inelastic creep contribution is of the standard power- lav form 

1/2 

with s^j as the stress deviator and = (3s^^Sj^^/2) ' as the 
effective stress. This gives a creep- rate uniaxial tension at 
a given temperature as € = B{T)<r'^, where B is a strong function of 
temperature. The coefficient of thermal expansion a is also a 
known function of temperature. The creep term in (2) 
characterizes secondary creep so that a primary creep contribution 
is not Included in (2). In principle there would be no difficulty 
In including more accurate creep descriptions; however, available 
data for sillco:, does not appear to justify a more complicated 
description at this point. The data does permit specification of 
B(T) and a(T). 

In addition to the above equations, the equations of In-plane 
equilibrium, ^ * C, must be specified. The equivalent 
variational equation of equilibrium is: 
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A 


snd no boundary terns are involved since the tractions vanish 
everywhere on the boundary. Given the stress- rate field, the 
stress field is related to it by the steady-state relation 

®lj " ^ .X 

Since the stress is zero at the melt interface, (5) implies; 

X 


(5) 


a ( 6 ) 

o 

This completes the set of governing equations. 

The numerical method uses a velocity- based finite element 
method to discretize the fields. The problem is inherently 
nonlinear and therefore an iteration procedure is employed to 
obtain the solution. In principle, the iteration scheme works as 
follows. 



1. At the end of any iteration an estimate of the stresses 
is available. This estimate is used to calculate an estimate of 
the creep strain-rates in each element for the next iteration. 

2. The finite element method is used to obtain the next 
improved estimate of the velocity field. The method employs a 
variational equation where the creep- rate terms and the thermal 
strains both enter as body force terms. 

3. Given the new estimate of the velocities, the new 
estimate of the total strain-rate is available and this is used in 
(2) to obtain a new estimate of the stress-rate. 

4. Given the new estimate of the stress-rate, the new 
estimate of the stresses are obtained from (6). The steps in a 
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'filvftn it* rat ion are now oonplete, and the next Iteration can atart 
over with step 1. 

Iteration is continued until convergence is obtained. The 
above scheme is modified in several nonessential ways to improve 
efficiency and stability. 

During the quarterly period starting from July 9, 1982, the 
existing computer code from the work in (1) and (2) was modified 
to be applicable to the ribbon growth problem. Efforts to check 
out the program on some trial elastic problems were started. 

These efforts will be continued into the next quarterly period, 
when we will also initiate calculations using realistic 
temperature distributions with creep deformations included. 

B. field Deterrojjsaljoa iD MfiviBS Sili£2D 

(R.O. otfll and J.P. Kalejs, Mobil Solar) 

In order to calculate the stress distribution and how it 

varies with growth speed in a solid, it is necessary to have a 

rather precise knowledge of the temperature distribution as a 

function of position and time. One possibility is that of 

measuring the temperature using thermocouples attached to the 

ribbon, but the large number of points required makes this 

experimentally impractical. Thus, we are left with the necessity 

of determining the temperature distribution by calculation and the 

use of experiment to validate the results by measuring a few 

selected points. 

The calculations have been made in two steps. The first was 
to use the surrounding temperatures and geometry to calculate the 
heat flux arriving at the surface of the ribbon. This was done by 


6 



OF POOR QUAumr 


••tUBlny th« furrounding turfae«» w«r« bI«ok bodl«i and 
integrating view factora to calculate the fliuc. 

The second step was to calculate the interaction of this flux 
with the silicon ribbon making use of the physical properties of 
'he material. Silicon is not a black body, but rather a gray body 
that is semi> transparent for energies below its band gap. Thus, 
it is necessary to calculate all energy interchanges as a function 
of wavelength and integrate to obtain the total net radiation 
flux. 


Some of the more important assumptions are the following. 


1. For photon energies less than the band gap. silicon 
is semi-transparent, but for photon energies greater 
it is opaque. 

2. Below the energy gap. the absorptivity and 
emisslvity are determined by the free carrier 
absorption and are functions of wavelength. Since 
at high temperatures silicon is intrinsic, the 
carrier concentration depends only on temperature. 

3. The emisslvity above the energy gap is one minus the 
reflectivity. The reflectivity at all wavelengths 
is a function of tne index of refraction and thus of 
wavelength. 

4. The energy gap and thermal conductivity are 
functions of temperature. 


5. The differential equation governing heat transport 
is that for a "thin” two-dimensional solid and is 
given by 


/ir ^ d /fc dT\ ^ - dT 

dr \ d 2 / dy \ dy/ di 


1 F(T,z.y) 
t 


where k is the thermal conductivity, v is the 
velocity in the z direction, C is the specific 
heat, t is the ribbon thickness and F(T,z,y), which 
is an integral over wavelength, is the difference in 
flux arriving and leaving the surface. 
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to oay. not nuoh can bo dont onolytioolly with ouch 
o highly non-Ilntoc oquotion ond, unless properly chosen » even 
nunericel techniques cen be unset isfectory. We heve adapted a 
technique called "quasi-linearization** and have obtained 
convergence in a relatively few iterations/ ' This procedure is 
closely related to the Hewton-Raphson method. 

In order to produce a condition where experiment and theory 
could be compared, a highly idealized situation was constructed in 
ribbon growth Furnace 17. A long slot (large enough to 
accommodate a S cm wide ribbon about 25 cm in length) wat milled 
in a block of fiberform insulation. Thermocouples were embedded 
in the walls to measure the surface temperature and two more 
thermocouples were attached to the ribbon. A plate at the bottom 
of the slot was heated to near the melting point of silicon and 
its temperature was monitored also. With the ribbon in the slot, 
and resting against the plate, thermocouple readings were taken. 
This was done for two different slot aspect ratios. Initial 
comparisons between calculations (in this case the ribbon velocity 
is zero) and experiment seemed to indicate that in the experiment 
the ribbon was 65^C too hot. This discrepancy was resolved when 
we determined that the thermocouples embedded in the walls were 
not recording the actual wall surface temperature because they 
were about 1 mm from the surface. They were low by about the 
right amount to explain the disagreement. 

The second calculation has been to model the expected 
temperature near the center of a ribbon for a typical **cold shoe** 
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'gtOMtry. Plgut* 1 thovt th« 9 «OMtry AodtXtd. Dit top rodlatlon 
•hiold (not shown), cold shoo and l^lntar cooling plata 
toMparataras wara all obtainad axparinantally fron naaaurasianta 
using tharmocouplas anbaddad in tha raspactiva oonponants. Thay 
thus dafina a ribbon anvironmant with prasoribad radiating 
anclosura wall tamperaturas. Tha linaar cooling plata tauparatura 
profila is typical of that naasurad for powar laval sattings usad 
in growth of 10 cm wlda ribbon. Tha cold shea, which is at 4S0^C 
and 0.3 cm long, is located 0.2 cm abova tha die. Wa assume th»t 
all tamperaturas vary linearly with position batwaan tha positions 
at which thay wara measured. A ribbon of thickness of 300 pm 
centered in the 0.2S cm gap is considered. Tha unknown 
temperatures in the model represent the dia-radiation shield gap 
(1500^0, tha radiation shield-cold shoe gap ( ilOO^C-SOO^C ) , and 
the cold shoe- linear cooling plata gap (700^C-1000^C). Tha 
sensitivity of tha ribbon temperature profila to each of these 
temperature parameters, as well as to tha enclosure size (ribbon 
to wall spacing) will be investigated in tha future. 

Tha calculations vara made with tha boundary conditions on 
tha ribbon that at tha die tha temperature is that of malting 
silicon (1412^0 and at tha exit to the furnace dT/dz * 0. Wa 
assume the horizontal flux is constant across tha ribbon width and 
at tha edges dT/dy * 0. 

Figure 2 shows tha ribbon tamparatura as a function of 
distance from tha die. Also shown is the tamparatura of tha 
environment. As expected, the tamparatura raacrss a minimum 
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LINEAR 

COOLING 

PLATE 


COLD 

SHOE 


DIMENSION IN : em 



TEMPERATURE IN: C* 


I200*-HI00* 


1200 * 


1 000^1 190* 
700* 1000* 
450* 

IIOO^COO* 


Fig. 1. '*ColU shoe” geometry used to calculate the tenqperature distrilHition 
in ribbon. 
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a 300 pm thick ribbon 
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'(1053^0 slightly above the coXd shoe after v ich it increases to 
aaxinum (1188^0 in the afterheater. The gradient at the 
interface of 980^C/ca will support a growth rate of about 3 
ca/ain. This is lower than the deaonstrated speed capability of 
the 10 ca wide ribbon system, which is estimated to be about A to 
4.5 cm/fflin at this thickness ribbon. The sensitivity analysis 
planned for investigating the influence of the unknown 
temperatures parameterizing the system should be able to determine 
if this simplified model can predict the experimentally observed 
speed capability, while further refinement of the heat transfer 
model is also planned. 
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